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Isospin- violation in the nuclear force is a subject of current interest. Significant advances in 



Abstract 

■ We calculate at two-loop order in chiral perturbation theory the electromagnetic cor- 

rections to the dominant two-pion exchange nucleon-nucleon interaction that is generated 
by the isoscalar irN contact-vertex proportional to the large low-energy constant C3. We 
find that the respective 2-7T7-exchange potential contains sizeable isospin-breaking compo- 
nents which amount to about —1% of the strongly attractive isoscalar central 27r-exchange 
potential. The typical value of these novel charge-independence and charge-symmetry 
breaking central potentials is 0.3 MeV at a nucleon distance of r = m^ 1 = 1.4 fm. Our 
analytical result for this presumably dominant 27T7-exchange interaction is in a form such 
that it can be easily implemented into phase-shift analyses and few-body calculations. 
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the understanding of nuclear isospin-violation have been made in the past years by employing 
methods of effective field theory (in particular chiral perturbation theory). Van Kolck et al. [1] 
were the first to calculate (in a manifestly gauge-invariant way) the complete leading-order 
pion-photon exchange nucleon-nucleon interaction. In addition, the charge-independence and 
charge-symmetry breaking effects arising from the pion mass difference 777,^-1- — m n o = 4.59 MeV 
and the nucleon mass difference M n — M p = 1.29 MeV on the (leading order) two-pion exchange 
NN-potential have been worked out in Refs. [2, 3]. Recently, Epelbaum et al. [4] have continued 
this line of approach by deriving the subleading isospin-breaking 27r-exchange NN-potentials 
and classifying the relevant isospin-breaking four-nucleon contact terms. Moreover, some next- 
to-leading order corrections to the ^-exchange nucleon-nucleon potential (those proportional 
to the large isovector magnetic moment k v = 4.7 of the nucleon) as well as effects from virtual 
A-isobar excitation on the ^-exchange interaction have been calculated recently in Ref. [5]. 

The long-range (pion-induced) isospin-breaking potentials found so far turned out to be 
rather weak. Typically, their values at a distance of r = m" 1 = 1.4 fm lie below 50keV in 
magnitude (see e.g. Tables I and II in Ref. [5]). More pronounced and therefore numerically 
significant are the isospin-breaking effects in the subleading 27r-exchange interaction (see here- 
fore Figs. 8 and 10 in Ref. [4]). As in the isospin-conserving case this feature can be traced back 
to the large magnitude of the low-energy constants C2, C3 and C4 which enter as parameters in the 
second-order chiral pion-nucleon Lagrangian C\^. According to the numerical investigations 
made in Ref. [1] , the inclusion of the (leading order) ^-exchange potential had negligible effects 
on the 1 Sq low-energy parameters and it led to only a tiny improvement in the fits of the NN- 
scattering data. With the partial cancellations through the next-to-leading order corrections 
[5] taken into account one may then conclude that the long-range ^-exchange interaction is 
too weak to explain the isospin-breaking NN-observables. One reason for this feature lies in the 
relative weakness of the l7r-exchange interaction itself. For example, its dominant (isovector) 
tensor potential amounts to only 9.2 MeV at a distance of r = m™ 1 = 1.4 fm. In comparison 
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to this the isoscalar central NN-attraction generated by chiral two-pion exchange comes out 
significantly larger (see e.g. Fig. 8 in Ref. [6]). Therefore one may expect that electromagnetic 
corrections to this strongest long-range NN-interaction, i.e. the chiral 27T7-exchange interaction, 
will give rise to more substantial isospin-breaking effects. It is the purpose of the present short 
paper to evaluate a selected class of (presumably dominant) two-loop 27T7-exchange diagrams. 
Now, since the isoscalar central NN-attraction comes mainly from the ttttNN contact-vertex 
proportional to the large low-energy constant c 3 , one naturally chooses for a first exploratory 
study the class of two-loop 27T7-exchange diagrams with exactly one such c 3 -vertex. The es- 
sential long-distance information about the two-loop diagrams is contained in their spectral 
function, which we will evaluate here analytically. 
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Fig. 1: Electromagnetic corrections to the dominant isoscalar central 2n -exchange NN-interaction 
generated by the 7nrNN contact-vertex —2ic^f~ 2 S ab q a ■ Diagrams with the contact-vertex at 
the right nucleon line and diagrams turned upside-down are not shown. The spectral function 
ImT{iji) is calculated by cutting the intermediate mvy three-particle state. 
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As a result, we do indeed find sizeable charge-independence and charge-symmetry breaking 
central potentials of magnitude 0.3 MeV at r = m~ l = 1.4 fm (for c 3 = —3.3 GeV -1 ). These 
seem to be the largest long-range isospin-breaking NN-potentials obtained so far. In order 
to test their phenomenological relevance they should be implemented into future phase-shift 
analyses and few-body calculations. As a result, we do indeed find sizeable charge- independence 
and charge-symmetry breaking central potentials of magnitude 0.3 MeV at r = m -1 = 1.4 fm 
(for c 3 = —3.3 GeV -1 ). These seem to be the largest long-range isospin-breaking NN-potentials 
obtained so far. In order to test their phenomenological relevance they should be implemented 
into future phase-shift analyses and few-body calculations. 

Let us start with recalling the chiral 27r-exchange NN-interaction. Its dominant contribution 
to the isoscalar central channel comes from a triangle diagram involving one 7nrNN contact- 
vertex with momentum-dependent coupling: —2ic 3 f~ 2 5 ab q a ■ qb- The corresponding potential 
in coordinate space (obtained by a Laplace transformation) reads [6]: 

Vt\r) = V (6 + UZ + 10Z2 + ^ + ' (1) 

where z = m % r. The occurring parameters are: qa = 1-3 (the nucleon axial-vector coupling 
constant), f n = 92.4 MeV (the pion decay constant), m n = 139.57 MeV (the charged pion 
mass) and the low-energy constant c 3 = —3.3 GeV -1 . The latter value is the average of c 3 = 
—3.2 GeV -1 and c 3 = —3.4 GeV -1 obtained by Entem et al. [7] and Epelbaum et al. [8, 9], 
respectively, in fits to NN-scattering phase shifts at next-to-next-to-next-to-leading order in 
the chiral expansion. The numerical values in the first row of Table I display the magnitude 
and r-dependence of this dominant isoscalar central 27r-exchange potential Vq\t). 

Now we will add to the 27r-exchange triangle diagram a photon line which runs from one 
side to the other. There are five positions for the photon to start on the left hand side and 
seven positions to arrive at the right hand side. Leaving out those four diagrams which vanish 
in Feynman gauge (with photon propagator proportional to g^ v ) we get the 16 representative 
diagrams shown in Fig. 1. Except for diagram (I) these are to be understood as being duplicated 
by horizontally reflected partners. A further doubling of the number of diagrams comes from 
interchanging the role of both nucleons. 

We are interested here in the coordinate-space potentials generated by the 27T7-exchange 
diagrams shown in Fig. 1. For that purpose it is sufficient to know their spectral function or 
imaginary part lmT(ifi). Making use of (perturbative) unitarity in the form of the Cutkosky 
cutting rule we can calculate the two-loop spectral functions as integrals of the (subthreshold) 
NN — > 7T7T7 — > NN transition amplitudes over the Lorentz-invariant 27r7 three-particle phase 
space [10]. In the (conveniently chosen) center-of-mass frame this leads to two angular inte- 
grations and two integrals over the pion energies. Due to the heavy nucleon limit (M N — > oo) 
and the masslessness of the photon (m 7 = 0) several simplifications occur and therefore most 
of these integrations can actually be performed in closed analytical form. For a concise presen- 
tation of our results it is also advantageous to scale out all common (dimensionful) parameters 
from the spectral function: 

t rpi ■ \ ucz9A m l Q ( V \ fo \ 
ImTM = ^W W' (2) 
and to work with the dimensionless variable u = f^/m n , where \x > 2m n denotes the 7m / ~f- 
invariant mass and a = 1/137.036 is the fine-structure constant. Without going into further 
technical details, we enumerate now the contributions of the 16 representative 27T7-exchange 
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diagrams shown in Fig. 1 to the dimensionless spectral function S(u). We find for u > 

S(u)V = (2 - r 3 - rl){% - 5m + 4 + ^- + - 1u(m - 1)1, 

12 u — 1 u J 

where r 3 2 denotes the third components of the usual isospin operators. 

,S(m) (ii) = 36m - 2m 3 - 48 — ln(u - 1) , 

u — 1 u 

S(u) m = 9m 3 - 12m 2 - 54m + 84 + -(14m 2 - 11 - 3m 4 ) \n(u - 1) , 

XL 

S(u) {w) = 7m 3 - 12m 2 - 26m + 44 + -(8m 2 - 5 - 3m 4 ) ln(w - 1) , 

xl 

16 

S(m) (v) = 16m 2 - 12m 3 + 56m - 80 (m 2 - l) 2 ln(u - 1) 

M 

4 44 16 2 f u / 2 dx u(x + y)- 



1 S(m)( vi )=3m 3 -8m 2 -22m + 48 + ^ + -1u(m-1) + -(m 2 -2) f —In 

M — 1 M M Ji y 

with the abbreviation y = \J x 2 — 1. 



l y u(x- y) 



»S(m) (vii) = 20M 2 -9M 3 + 22M-52 + -(5M 4 -16M 2 + ll)ln(M-l) 

M 

8. 9 ^ W2 (mi — 1)(4 — m 2 — lux) u(x + y)-l 

+ — (u 2 -2)<b dx- p -In-) , 

Ji (m — 2x)y u(x — y) — 1 



u 



S(uY Ym) = (2-r 3 -r 3 )|3M 2 + 6M-ll+ (m 3 -2m-^1u(m-1) 

- 3 4 - J- + 1 (m 2 - 2) rt, il^+^-n 

2 m — 1 m Ji y u(x — y) — 1 ) 

S(u)W = (2 + r 3 + r 3 ){^-5M + 4+^- + -ln(M-l) 

K Z: XL X XL 

+1(26 - 5m 2 )v/m^4 - 8 ln M + v/ ^n , 
3 V ; m 2 J 

,S(m)( x ) = (2 + r 3 + r 3 )|M 2 -^ + y-5+(M 3 -4M+^ln(M-l) 

8,, 2 , r- / 3 in 8\ m + y/u 2 - 4 

+ -(1 - m ) v m - 4 + 4m 3 - 10m + - In — 
3 V m/ 



+-(u 2 -2)f — ux/ + (1 - w 2 ) In y - 

u Ji u — 2x1 u — x — y \ } 
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S(u)^ = S'( w )( x ) + (2 + 7f + r 2 3 ){« 2 -« 3 +(6«-« 3 -^)lii(«-l) 

1 2 . 2 . f u / 2 1 ux — 1 m(x + m) - 1 "I 

+3m-1 - + -m 2 -2 / rfx In-) ^ , 

m — 1 m Ji y u(x — ?/) — 1 J 



S(u)^ = ( T ? + Tl + 2T?Tl)l--u 2 - — + 5+Uu-u 3 --)\n(u-l) 

14 2 V u J 

-3u 2 V^4 + Uu 3 -6u+ 8 -) In U _±^E± 
V uj 2 



+-«-2 — ™/+l-u 2 ln ^ , 14 

u Ji u — 2x L -u — x — yJJ 

,S( M )( xni ) = S( M )( xn ) + (r 3 + r| + 2r 1 3 r 2 3 ) |2m 2 - y + 3« - 10 + ^2 M 3 - 8u + ^ ln( M - 1) } , ( 15) 

S(u)^ = (3 + r 3 + r\ - r 3 r|){ - 3 ) v^4 + 1 In " + V f^ | , (16) 

%)( XV » + %)( XVI > = 0. (17) 

The exact cancellation between diagram (XV) and the irreducible part of diagram (XVI) has 
the same kinematical reason as the cancellation between the crossed box and (the irreducible 
part of the) planar box diagram, which has been discussed in detail in Ref. [6]. We have 
carefully checked gauge-invariance. The total (dimensionless) spectral function S(u) stays 
^-independent when adding a longitudinal part to the photon propagator: g^ u — > g^ u + 
^k^k v . This property holds already separately in the subclasses of diagrams: (I) + (VIII), 
(II) + (IV) + (VI), (III) + (V) + (VII), (IX) + (X) + (XI), (XII) + (XIII) and (XIV). The "encircled" 
integrals appearing in Eqs. (9,12,14) involve the following regularization prescription: 

u/2 dx fa) f u/2 ^ f(%) - /(«/2) 



i u — 2x 



^tw-twv (18) 

J 1 XL ZiJj 



This regularization prescription eliminates from some diagrammatic contributions ((V), (VII), 
(X), (XI), (XII), (XIII)) to the spectral function an infrared singularity arising from the emission 
of soft photons (NN — > 7T7r7 so f t ). In the 27r7 phase space (an oval-shaped region in the cjicj 2 - 
plane spanned by the pion energies) the singularity is located at that extremal boundary point 
{uj\ t 2 = /x/2) where the photon energy /j, — U\ — 002 becomes identical to zero. The singular 
factor (u — 2x) _1 in Eq.(18) stems from a pion propagator. The regularization prescription 
defined in Eq.(18) is equivalent to the familiar "plus" -prescription [11] employed commonly 
for parton splitting functions in order to eliminate there an analogous infrared singularity due 
to soft gluon radiation. 1 We note as an aside that the non-elementary integrals (J^ 2 dx . . .) 
in Eqs. (8-14) can be solved in terms of dilogarithms and squared logarithms of the argument 
(u + Vu 2 - 4)/2 [12]. 

Let us first discuss some generic properties. All contributions to the dimensionless spectral 
function S(u) vanish at the 7T7T7-threshold u — 2. This is markedly different from the original 
one-loop two-pion exchange triangle diagrams which possess the discontinuous spectral function 
S(u)^ 2 ^ = —12n{au)~ 1 {u 2 — 2) 2 9{u — 2) with a non-vanishing threshold value. The leading 
threshold behavior of the total dimensionless spectral function S(u) is: S(u) = 64(1 + r 3 )(l + 
rf) (1— In 2) \/u — 2 + 0{u — 2) and it comes exclusively from the four diagrams (X)-(XIII). Such 
a square-root like growth implies an asymptotic tail of the coordinate-space potential of the form 
e -2m w r r -5/2 f Qr r _^ ^ whereas at short distances r — > the potential develops a r~ 6 In 2 (m n r) 

1 It would of course be desirable to extend the present calculational framework (by including various radiation 
diagrams etc.) such that the overall infrared finitcness could be demonstrated in detail. This goes beyond the 
scope of the present work and therefore we stay with the physically well-founded regularization prescription 
Eq.(18). 
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singularity, somewhat stronger than that of Vq (r) written in Eq.(l). At any intermediate 
nucleon distance r the 27T7-exchange central potential in coordinate-space V^ 2nJ \r) can be 
calculated from the spectral function lmT(ifi) via a modified Laplace transformation: 

Vt"\r) = --Vf^^ImT( V ) 

Z7l T J2m- J! 

= n°\r)+T?T!V^\r) + (T? + Tl)v£ sh \r), (19) 

where we have given in the second line of Eq.(19) its decomposition into isospin-conserving 
(0), charge-independence breaking (cib), and charge-symmetry breaking (csb) parts. We note 
again as an aside that the (partial) contributions to the potential V^ 27r7 ^(r) coming from the 
terms in Eqs.(3-16) without integrals can be represented by the exponential function e~ 2z , 
the exponential- integral function e~ z E 1 (z), and two modified Bessel functions K 0jl (2z), each 
multiplied by a pure polynomial in \jz of degree 6 or lower, where z = m w r [12]. 

We are now in the position to present numerical results for the dominant 27T7-exchange 
nucleon-nucleon potentials. The numbers in the second, third, and fourth row of Table I show 
the dropping of the central potentials V^°\r), V^ ah \r) and VQ Sh \r) with the nucleon distance 
r in the range 1.0 fm < r < 2.1 fm. Note that all potentials are given in units of MeV. 
One observes that the repulsive charge-independence breaking (cib) potential and the repulsive 
charge-symmetry breaking (csb) potential are approximately equal and about half as strong 
as the isospin-conserving one. Their magnitudes (for example, 0.3 MeV at a distance of r = 
1.4 fm) are exceptionally large in comparison to all so far known isospin-breaking NN-potentials 
generated by long-range pion-exchange. The dominant contributions to VQ lb \r) and Vc° sb \r) 
come from the four diagrams (X)-(XIII). This selective feature is consistent with the above- 
mentioned leading threshold behavior of the spectral function S(u) and the associated large-r 
asymptotics. 



r [fm] 


1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2.0 2.1 


v c 

V C 
T/ (cib) 

T/ (csb) 
V C 


-221 -121 -69.3 -41.4 -25.6 -16.3 -10.6 -7.05 -4.78 -3.30 -2.31 -1.64 
6.47 3.25 1.72 0.953 0.547 0.325 0.198 0.124 0.079 0.051 0.034 0.023 
3.54 1.78 0.948 0.526 0.303 0.180 0.110 0.069 0.044 0.029 0.019 0.013 
3.57 1.80 0.956 0.530 0.305 0.182 0.111 0.070 0.045 0.029 0.019 0.013 



Table I: The dominant isoscalar central 2n -exchange NN-potential Vq (r), and electromagnetic 
corrections to it, as a function of the nucleon distance r. The values in the third and fourth row 
correspond to the is o spin-violating potentials V c ah \r) and V c csb \r). All potentials are given in 
units of Me V. 

It is also instructive to compare the isospin-violating 27T7-exchange potentials with the 
isospin-conserving 27r-exchange potential V c 2n \r). Their relative ratio of about —1 percent 
is consistent with the usual rule of thumb estimate, namely a/n ~ 1/430 times a numerical 
factor. In the present case this numerical factor is not just 1 but of the order 5, due to the 
large number of contributing diagrams (in total there were 54 non- vanishing diagrams). 

We have also computed the 27T7-exchange central potentials proportional to the other large 
low-energy constant c 2 = 3.3 GeV -1 . At a reference distance of r = 1.4 fm one finds the values 
V c °\r ) = -8.54 keV, V c cih \r ) = -23.4 keV and % csh \r ) = -14.8 keV. These potentials 
are in fact more than an order of magnitude smaller than their counterparts generated by the 
C3-vertex. Furthermore, as expected the 27T7-exchange central potentials arising from the c\- 
vertex are small. With Cl = -0.8 GeV" 1 we get V^\r ) = 23.4 keV, V^ db) (r ) = 8.40 keV and 
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Vq (r ) = 9.78 keV. This confirms our assumption made about the dominant 27T7-exchange 
NN-potential. 

In summary we have calculated in this work the electromagnetic corrections to the dominant 
two-pion exchange nucleon-nucleon interaction that is generated by the isoscalar ttN contact- 
vertex proportional to the large low-energy constant c 3 = —3.3 GeV -1 . We have evaluated ana- 
lytically the spectral function of a large class of 70 two-loop 27T7-exchange diagrams. The result 
for the total spectral function is manifestly gauge-invariant. An infrared singularity related to 
soft photon emission which occurs in some contributions has been regularized by a prescription 
analogous to the one employed commonly for parton splitting functions. As a major result 
we have found that the 27T7-exchange potential contains sizeable isospin-breaking components 
which amount to about —1% of the strongly attractive isoscalar central 27r-exchange potential. 
A typical value of the charge-independence and charge-symmetry breaking central potential is 
0.3 MeV at a nucleon distance of r = m" 1 = 1.4 fm. Our analytical result for this novel and 
exceptionally large isospin- violating long-range nuclear force has been presented in a form such 
that it can be easily implemented into phase-shift analyses and few-body calculations. For 
completeness one should also calculate the 27T7-exchange diagrams proportional to the remain- 
ing large low-energy constant c 4 ~ 3.4 GeV -1 which generate (isospin- violating) spin-spin and 
tensor potentials. Some of the technical advantages arising from the Lorentz-scalar nature of 
the c 3 -vertex are then, however, no more available. The issue of (overall) infrared finiteness 
should also be further explored. Work along this line is in progress. 
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